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1. Introduction 

The development of large momentum transfer (LMT) mirrors and beamsplitters for 
atoms is a promising path for increasing precision in future sensors based on atom 
interferometry PQ El E]- In particular, Bragg diffraction of atoms from a light pulse 
seems to be a strong candidate for LMT atom-optical elements, and has received much 
theoretical [U O |6l [7] and experimental attention [HI [9l [TOl HH [12]. However, much 
theoretical work on Bragg diffraction assumes a plane-wave atomic source. In this 
paper, we relax this assumption, and theoretically investigate how the momentum width 
of an atomic source affects the efficiency of Bragg beamsplitters and mirrors and, more 
generally, the signal-to-noise ratio (SNR) of a Mach-Zehnder (MZ) interferometer. In 
particular, we demonstrate that the momentum width places fundamental restrictions 
on the Bragg mirror efficiency, and consequently limitations on the interferometer 
sensitivity. In addition, limitations on laser power and spontaneous loss place further 
restrictions on efficiency and sensitivity that can only be overcome by severely narrowing 
the momentum width. The analysis presented in this paper gives one the means 
of optimizing Bragg mirror and beamsplitter pulses, and furthermore allows one to 
determine the source type, momentum width and momentum transfer that maximize 
the signal-to-noise ratio of a Bragg pulse atom interferometer. 

To date, cold atom inertial sensing has been dominated by thermal sources typically 
at a temperature on the order of 1 fiK or less. Whilst such sources can be created with a 
high atom flux, they typically have a large momentum width on the order of 1 hk, where 
k = |k| is the wavenumber of the laser light. In experimental work on LMT mirrors and 
beamsplitters carried out by a number of groups, it has been necessary to velocity select 
the source cloud (reducing the atom flux) prior to beamsplitting in order to realize high 
fidelity beamsplitting and reflection [U [TTl [13]. This demonstrates that, practically, a 
narrower momentum width is preferable even though it leads to a reduction in total 
atom flux. 

There has been some interest in understanding whether it is advantageous to use 
a narrow momentum source, such as a Bose-Einstein condensate (BEC), for atom 
interferometry. For example, Kozuma et al. [H] and Torii et al. [13] have used low order 
Bragg diffraction to transfer a condensate to a target momentum state with an efficiency 
close to 100%. Hughes et al. [lOJ considered how the fidelity of Bragg beamsplitters 
and mirrors depends upon the atomic velocity of a BEC, and hence showed that the 
narrow momentum distribution of a BEC allowed for very efficient 1st, 2nd and 3rd 
order Bragg beamsplitters and mirrors. More recently, the phase sensitivity of the atom 
interferometer developed by Impens et al. [IHl |T6] , which involves the levitation of a 
condensate, has been shown to improve when the atomic source is progressively narrowed 
through filtering with multiple wave atomic interferences. Debs et al. [12] constructed 
a BEC gravimeter, utilizing Bragg LMT mirrors and beamsplitters, that had superior 
fringe contrast to a similar thermal source gravimeter. 

Whilst such results illustrate the importance of a narrow momentum width source. 
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one must recall that current Bose-condensed sources have a lower atom flux than state-of- 
the-art thermal sources. Indeed, a state-of-the-art BEC machine that produces 2.4 x 10^ 
atoms/s [IT] still has a flux 25 times lower than the best published state- and velocity- 
selected thermal cloud used in a gravimeter ^18]. It has thus been unclear whether the 
advantages of a narrower momentum width compensate for the reduced atom flux. We 
quantitatively address this question, and accommodate both flux and momentum width 
considerations by optimizing the only parameter that matters: the signal-to-noise ratio. 

Previous experimental and theoretical work has also suggested that as the order 
of the Bragg process is increased (which increases the momentum imparted during the 
process), a narrower momentum width is required in order to realize high efficiency 
beamsplitting and mirror processes. For example Catallioti et al. [19] reported 
experimental and theoretical results indicating that the linewidth of multiphoton 
transitions decreases as the reciprocal of the order of the process, at least in the 
perturbative regime where the excited state population of the two-level system is kept 
small. For a condensate, Wu et al. [20] showed that the transfer efficiency of a particular 
two-pulse Bragg beamsplitter decreased as the Bragg order increased. Interestingly, 
our analysis shows that the maximum possible transfer efficiency of Bragg mirrors and 
beamsplitters does not degrade at higher Bragg orders for clouds of momentum width 
less than approximately 0.25 hk, on the proviso that one is not limited by laser power. 
Whilst the maximum transfer efficiences of larger momentum width clouds do vary 
with Bragg order, this dependence is small and is almost negligible for n > 5. However, 
when the two-photon Rabi frequency is restricted our analysis reveals that the efficiency 
of higher order Bragg processes is reduced by a factor that strongly depends upon 
the momentum width of the cloud. As the Bragg order increases, this degradation 
in efficiency becomes more severe and can only be overcome by increasingly narrow 
momentum widths and longer pulse durations. This demonstrates that the practical 
implementation of high fidelity LMT Bragg mirrors and beamsplitters necessarily 
requires a narrow momentum width cloud and a large laser intensity. 

Our paper is structured as follows. Section [2] briefly introduces a MZ atom 
interferometer and motivates the interest in LMT for atom interferometry. Section [3] 
introduces a model that describes Bragg diffraction, and develops a criterion that can be 
used to optimize Bragg mirror pulses. In section 4.1 , numerical simulations using this 
model and criterion demonstrate that: (a) the momentum width places a fundamental 
limit on the reflection efficiency of a Bragg mirror; (b) for unconstrained two-photon 
Rabi frequency, the maximum efficiency of a mirror pulse does not markedly change 
with increasing Bragg order; and (c) the optimal two-photon Rabi frequency scales 
roughly quadratically with the order of the Bragg process. Qualitative features of this 



numerical analysis are explained using a simple two-level model in section 4.2 In 
section 4.3, it is shown that when the two-photon Rabi frequency is limited (due 



to limitations on laser power and spontaneous loss), (a) the efficiency of a Bragg 
mirror drastically decreases above some Bragg order; and (b) achieving the same 
mirror efficiency at a higher Bragg order requires a reduction in the momentum width. 
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The latter half of the paper is concerned with performing a similar analysis for a 
MZ interferometer constructed with Bragg mirrors and beamsplitters. Specifically, 
section 5.1 constructs an efficiency criterion for a MZ interferometer, thus allowing one 
to optimize both Bragg beamsplitters and mirrors, section 5.2 contains a qualitative 
analysis of a MZ interferometer with the two-level model, and section 5.3 presents the 
quantitative results of the numerical optimization, with and without constraints on the 
two-photon Rabi frequency. The results mimic that of the Bragg mirror, showing that 
the momentum width fundamentally limits the maximum SNR for a phase measurement 
with a MZ interferometer, and particularly so when the two-photon Rabi frequency is 
constrained. In addition, a constraint on the two-photon Rabi frequency implies that, 
above some Bragg order, the SNR actually worsens. Finally, in section 5.4 the optimal 
Bragg order and (for a thermal source) velocity selection that maximizes the SNR of a 
MZ interferometer are calculated. A comparison of current state-of-the-art thermal and 
Bose-condensed sources demonstrates that, in principle, current Bose-condensed sources 
are competitive with state-of-the-art thermal sources. It is argued that future increases 
in atomic flux and the practical advantages of a narrow momentum width source imply 
that Bose-condensed sources have a signiflcant role to play in the next generation of 
precision sensors based on Bragg pulse atom interferometry. 



2. Background 

A comprehensive review of atom interferometry is given in [21]. We consider an atom 
interferometer in the Mach-Zehnder (7r/2 — vr — 7r/2) conflguration, where nth order 
Bragg transitions are used as beamsplitters (7r/2 pulse) and mirrors (vr pulse). Such an 
interferometer operates by coupling two momentum states \pq) and \pq + 2nhk), where 
k = |k| is the wavenumber of the laser light and n is the (integer) order of the Bragg 
process. If the acceleration of the atoms relative to the laser is uniform, then the free 
evolution phase in each arm is identical, and the only interferometric phase contribution 
is from the atom-light interaction. In particular, it is straightforward to show that the 
interferometric phase shift is proportional to the acceleration [22J: 

$ = n(0i - 202 + 03) = 2nk ■ aT^, (1) 

where 0j is the optical phase of the ith Bragg pulse at the point where the atom interacts 
with the pulse, relative to some reference position, a is the uniform acceleration (due 
to gravity, say) relative to that reference position, and T is the time between pulses. 
If we scan the phase $ and measure the relative population in the momentum state 
\po + 2nhk), then we obtain a signal proportional to (1 + C cos $) /2, where A^ is the 
atomic population in the two momentum states and C is the contrast of the fringes. 
Therefore, for operation at mid-fringe and small phase shift the SNR at the shot 
noise limit is 

SNR oc Cy/N {2nhkT^) ■ a. (2) 
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(a) Ideal Bragg mirror process (plane-wave input, effective two-level system) 



<1 00% transfer 




(b) Realistic mirror process (general initial momentum distribution, multi-level system) 



Figure 1. Energy iJ^jn vs momentum p diagrams illustrating the difference between 
an ideal Bragg mirror process and a realistic diffraction process. Here = Huir is the 
recoil energy. Without loss of generality, these diagrams show the case where po — 0. 
An ideal Bragg process is shown in (a). In the realistic process (b), the initial atomic 
cloud has a finite momentum width. In this case it cannot be transferred to the target 
state with 100% efficiency, and furthermore there is loss to diffraction orders other 
than the two interferometer modes. 

The motivation for LMT with Bragg transitions is now clear: the signal can be increased 
in direct proportion to increasing the momentum transferred between the two states (i.e. 
increasing the Bragg order n). 

In an ideal Bragg process, an atom in a single momentum state \po) is coupled to a 
final momentum state |po + 2nhk). Specifically, for nth order diffraction, n photons at 
frequency ui and wave vector k are absorbed and stimulated to emit into frequency uj2 
with wavevector — k, leading to a transfer of momentum 2nhk and an energy of nhAun 
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[cf. figure l(a)| . However, as shown in figure 1(b), in a realistic system the atomic 



cloud will initially be in some momentum distribution of width a. Hence we cannot 
expect that, in general, perfect vr and 7r/2 pulses exist. Furthermore, realistic diffraction 
processes couple the initial momentum state \po) to a range of momentum states. This 
results in a loss of flux from the interferometer, which we would like to minimize. Indeed, 
as shown by Miiller et al. [6], practical limitations on pulse duration require operation in 
the multi-level quasi-Bragg regime, which is the intermediate regime between the long 
interaction Bragg regime and the short interaction Raman-Nath regime. The theoretical 
model introduced in the next section places no restriction on the momentum width of 
the atomic cloud, nor the number of momentum states that take part in the diffraction 
process. 

3. Theoretical model 

We consider a semiclassical model of Bragg scattering, where a cloud of two-level atoms 
interacts with a classical light field. This model has been successfully applied to similar 
Bragg scattering problems [23], and is therefore suitable for our purposes. Whilst figure |2] 
shows Bragg diffraction in the frame where a stationary atom interacts with a moving 
optical lattice, it is more convenient to work in the frame where the lattice is stationary. 
In this case, an nth order Bragg resonance corresponds to a resonant atomic momentum 
of —nhk, for integer n. At this point, we assume an initial momentum state of {—n+6)hk. 
Here 6hk can be thought of as the momentum detuning from resonant momentum. 
Provided the laser is far detuned from the excited state, the atomic excited state can 
be adiabatically eliminated as shown in [6]. In this limit, the dynamics are completely 
described by the Schrodinger equation for the ground state wavefunction g[z,6,t): 

thg{z,S,t) = -^ ^'^g/'^^ + 2^^(t) cos'ikz)giz,S,t), (3) 

where M is the mass of the atom, z is the position co-ordinate along the direction of 
the standing wave and fl(t) is the two-photon Rabi frequency. We have assumed that 
the atom-atom interaction energy is negligible and can be ignored. This is true for most 
atom interferometer experiments, including those using expanded BEC sources [121 El] , 
where these interactions would give an undesired systematic effect. 

Since the light potential is periodic, we apply Bloch's theorem and express g{z, S, t) 
in terms of eigenfunctions of constant quasimomentum: 

oo 

g{z,5,t)= C2™(5,t)e^(2'"+^)'=^ (4) 

m,=— oo 

Substituting Q into ^ gives the following set of ordinary differential equations 
(indexed by integer m): 

iCm{S, t) = Ur{m + 6fCm{S, t) + ^ [^](^)c„+2(5, t) + (t)c„_2 (5, t)] , (5) 

where Ur = hk'^/2M is the atomic recoil frequency. For the 780 nm D2 transition of 
s^Rb, ~ 27r X 3.8 kHz. 



Why momentum width matters for atom interferometry 



7 



In what follows, we will consider Gaussian pulses, as they give lower losses to other 
diffraction orders than square pulses [6]. For a single pulse, 

Q{t) = Qexp{-ty2T^) , (6) 

where fl is the amplitude of the pulse and r is the pulse duration. For a MZ 
interferometer, there are three Gaussian pulses, separated by time T, and the phase 
difference between the two arms is included by introducing a phase shift on the third 
pulse: 

n{t) = fibsexp {-ty2Tl) + n^exp (-(t + r)V2r2) 

+ fibs exp {-{t + 2T) V2r2, + tcj)) . (7) 

Here fibs (fim) and rbs (rm) are the pulse amplitude and duration for the beamsphtter 
(mirror) , respectively. 

3.1. Efficiency criterion for a mirror pulse 

Initially, we assume the atomic cloud to be in the state 

\i^i)= [ pmi-n + S)hk)d6, (8) 



where p{6) is the cloud's momentum distribution initially centred at momentum —nhk. 
We assume p to be real, and note that it must satisfy J p{S)^ d6 = 1 in order for \ipi) to 
be correctly normalised. 

For an nth order Bragg transition a mirror pulse ideally maps lipi) to the final state 

l^ideai) = j p{S)\in + 6)hk)d6. (9) 

That is, a distribution p of momentum states centred at momentum —nhk has been 
mapped to the same distribution p of momentum states centred at momentum nhk. 
However, in general we will get a final state of the form 

\^) = Y.j Cm{S,tf)\{m + 6)hk)dS, (10) 

m 

where the coefficients Cm{5,t) are determined by solving the set of equations ([s]), and 
if is some time well after the pulse has been turned on and off (i.e. fi(t/) ~ 0). We 
characterize the difference between the state we get and the state we want through the 
fidelity: 

= I (^Ideall^) p. (11) 

By direct substitution of ^ and (10) into (11 ), and application of the Cauchy-Schwartz 
inequality, we can prove that 



F^< J p{6y\cn{6,Tf)\US. (12) 

If we assume that the relative phase between different momentum components of the 
cloud is approximately zero, then this upper bound will approximately be the fidehty. 
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Thus, if we take p{6) to be a normalised Gaussian of width a, the efficiency criterion 
for a Bragg mirror is 

F.^J y=^"''^''^' \Cn{6,Tj)f d6. (13) 

Given that e"'^^/^'^^ \cn{S,Tf)f /y/2no^ is the population transferred to the state |(n + 
6)hk), we can see that F-j^ is just the total population that is transferred to a distribution 
centred at momentum nhk. Clearly, the best Bragg mirror will maximize transfer 
to this distribution, hence maximizing the fidelity. 

The efficiency criterion defines our optimization problem. Specifically, our goal is, 
for different Bragg orders n and momentum widths a, to find the optimal two-photon 



Rabi frequency flopt and pulse duration Topt that maximize criterion (13). 

4. Analysis for a Bragg mirror 

4.I. Numerical simulation with unconstrained power 

For different Bragg orders n and momentum widths a, the optimal pulse amplitude 



Qopt and duration Topt were found by maximizing the fidelity (13). was computed by 
numerically solving a truncated set of equations (|5| for various 5, and r. Momentum 
states larger (smaller) than mmax (— ''^max) were truncated, as they were not populated 
during the diffraction process. Given that the initial and target states have kinetic 
energy n?hjjjr and the potential energy due to the optical lattice is Wt, the highest energy 
state m to be populated is given by n'^ujr + ^ ~ m?ujr. Including the nearest two or three 
higher energy states to m, we thus set mmax to the integer part of {^JVl/ ujr + rfi + 6) if 
n was even, or to the integer part of {^JVt/ur + n? + 7) if n was odd. 

For simplicity, we restricted the optimization for r to be within the first Rabi cycle. 
This is not an unreasonable restriction, given that an interferometer requires pulse 
durations much shorter than the interrogation time (i.e. the time between successive 
pulses). 

Figure [2] displays plots of the maximized fidelity F^(fiopt, Topt) as a function of the 
momentum width a for a variety of Bragg orders n. The first point to note is that for 
a fixed n, monotonically decreases with increasing a. Thus the size of the cloud's 
momentum width places a fundamental limitation on the maximum transfer efficiency 
of the mirror. Secondly, for fixed a, increasing n does not necessarily change the 
maximum F^ obtainable. A clear example of this is at a = 0.1 hk, where the maximum 
transfer efficiency is roughly F^^ ~ 0.9 for n > 2. Furthermore, as shown in figure [2]^b), 
Fni^opt, Topt) is almost entirely independent of n for higher order Bragg processes. The 
caveat, however, is that for larger n, maximal requires a larger Q, and hence a larger 
laser intensity (or smaller detuning). Our simulations show that, irrespective of a, Qopt 
scales roughly quadratically with n (see figure [3]). That is, optimal refiection occurs 
when Q is on the order of the kinetic energy, n^w^. This is consistent with the scaling 
in the analytic limit of the two-level Bragg regime of ~ ((Sw^)"^^ [{n — 1)!]^) ^ [cf. 
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Figure 2. Plots of the maximum fidelity -^^(f^opt, ''opt) as a function of momentum 
width a. Different curves correspond to different order Bragg processes. The curves 
have been separated into two plots to better allow for discrimination between them. 
Note that the actual data is given by the points; the lines are simply to guide the eye. 



section 



4.2, and in particular (15)]. As an aside, the typical values for Topt ranged from 



0.08 - 0.5 w"^ (see figure |3|). 

As a simple application, we compared our calculations of the mirror efficiencies 
to measured values in the literature. In Miiller et al. [H] mirror efficiencies between 
80 — 90% were reported for Bragg orders n < 9 with a source of momentum width 
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Figure 3. Plots of (left) optimal two-photon Rabi frequencies and (right) optimal 
pulse durations as functions of Bragg order. These optimal pulse parameters give the 
maximum fidelity achievable for an nth order Bragg mirror process and a cloud of 
momentum width cr (cf. figure l2|. Lines joining the points are simply to guide the eye. 



a ~ 0.13 hk. This is consistent with our theory's prediction of ~ 0.85 — 0.9, and also 
suggests that Miiller et al. were operating close to the fundamental mirror efficiency 
limit. Similarly, Chiow et al. [25] observed 6 hk reflections of a BEC with efficiency 94%. 
From [2S], we estimate the momentum width of their expanded cloud to be a ~ 0.05 hk, 
which places an upper bound of approximately 95% on the mirror efficiency. In both 
cases, significant improvements to the mirror efficiency would only be possible by a 
further narrowing of the source momentum width. 

Another stark feature seen in figure [2]^a) is the 'cusps', where at a particular a 
the plot undergoes a sudden change in slope. This cXcusp corresponds to a change in 
the optimal operating regime. At a < (Xcusp the atomic population is almost entirely 
divided between momentum distributions around the initial and target states, | — nhk) 
and \nhk) respectively. For a > (Xcusp the optimal pulse results in non- negligible loss to 
other diffraction orders. 

There has been some interest in constructing LMT atom-optical elements with a 
series of lower order Bragg transitions, rather than a single higher order transition [25] . 
Our results suggest that, provided one has ample laser power, a single nth order Bragg 
pulse will be more efficient than an mth order pulse applied q times, where qm = n. 
For instance, for an atomic cloud with momentum width 0.2 hk an n = 10 Bragg mirror 
process has a maximum fidelity of approximately 0.77. Thus about 23% of the atoms will 
not be reflected, and in the simplest case would leave the interferometer. In contrast, a 
series of flve n = 2 pulses, where each pulse also has a maximum fldelity of 0.77, would 
have an effective maximum fldelity on the order of 0.77^ ~ 0.27. Ignoring reductions 
in contrast, this loss of atom flux would reduce the SNR by a factor of two. Similar 
losses in atom flux would affect multiple Bragg pulse beamsplitters. This degradation in 
mirror and beamsplitter fldelity could be somewhat ameliorated by using an extremely 
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Figure 4. Plot of the efficiency criterion as a function of the w parameter [cf. (17)] 
Note that the function monotonically decreases with increasing w. 



narrow source, such as an atom laser. Nevertheless, all other considerations being equal, 
a single Bragg pulse will give a more efficient 2nhk reflection than a series of lower order 
Bragg pulses. 

4-2. Two-level model for Bragg mirror 

We can qualitatively understand the dependence of the fidelity on the momentum width 
a, Bragg order n and two-photon Rabi frequency fl by considering the following simple 
model. Consider Bragg diffraction from a square pulse in the Bragg regime, where 
requirements of energy conservation are so severe that only two states are involved 
in the diffraction. More formally, the intermediate states between \{—n + 6)hk) and 
\{n + 5)hk) are adiabatically eliminated, and so the set of equations ^ reduce to [H] 

ic-n{5, t) = -2urn6cn{S, t) + -^Cn{5, t) (14a) 

iCn{S,t) = 2urn6cn{S,t) + ^c_n{S,t), (146) 

where Qes is the effective (or 2n-photon) Rabi frequency. In the regime a -C hk, fles is 
approximately independent of 6 and is given by 

"^^^=(8.% [(n-l)!f ^''^ 



Solving (14a) and (14&), with initial condition C-n{5,0) = 1 and Cn{5,0) = 0, for the 



target state population gives 



Comparing with the standard two-level atom, it is obvious that AuruS behaves precisely 
as a detuning from resonance. Given that most of the cloud is centred at 5 = 0, and 
we can get the maximum population transfer on resonance, we choose a pulse time 
t = n/Qca- 
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Substituting (16) into the efficiency criterion (13) gives 



where we have transformed to the dimensionless variables 
a; = — ; ^ = "PS • 

i ieff i 'eff 

We can now see that in this two-level regime the only important parameter governing 
the efficiency of a Bragg process is w. A plot of as a function of w is shown in 
figure |4| Noting that F-j^ monotonically decreases with w, it is clear that the two- level 
model of the Bragg mirror predicts the following: 

• For a fixed figg and n, an increase in momentum width a increases w, and thus 
decreases the efficiency of the process. 

• For a fixed a and Qes, an increase in the order of the Bragg process n gives lower 
efficiency mirrors. 

• If it is possible to increase figfr, then decreases in efficiency due to increases in n 
can be nullified. That is, for higher order Bragg processes and larger momentum 
width clouds, we can achieve a low w (and hence a high F^) by increasing the laser 
power. 

This simple model qualitatively captures key features of our comprehensive numerical 



simulation of the Bragg mirror (cf. section 4.1). However, it does not quantitatively give 



the correct maximum fidelities, nor the correct optimal pulse amplitudes and durations. 
Furthermore, the model predicts that we can compensate for increases in a by increases 
in f2efj. Figure [2] shows that this is manifestly false for a realistic Bragg reflection. 
These disagreements are not surprising, given that the model only has validity in the 
regime where (a) the momentum width is much narrower than hk; and (b) for each 6, the 
intermediate states between \ {—n + 6)hk) and \ {n + 6)hk) are far detuned from the Bragg 
resonance (this bounds the effective Rabi frequency to Qes ^ ScOr{n — l)"/[(n — 1)!]^). 

4-3. Optimization with limitations on the two-photon Rabi frequency 

The numerical results of figure [2] show that the maximum fidelity remains virtually 
unchanged as one increases the Bragg order, on the proviso that one increases the two- 
photon Rabi frequency such that it is on the order of f2 ~ n'^Ur- However, there will 
always be practical limitations on laser power and the loss that can be tolerated due 
to spontaneous emission. Both these constraints will place a clamp on the two-photon 



Rabi frequency. The two-level model from section |4.2| predicts that if there is a clamp 
on QcS, then increasing the Bragg order n past a certain point results in an increase 
in w, and hence a rapid decrease in the mirror fidelity. A comprehensive numerical 
analysis shows this result to be true in general, and furthermore shows that a clamp 
on the two-photon Rabi frequency limits high Bragg order, high fidelity reflections to 
narrow momentum width sources. 
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We can estimate the Bragg orders where the maximum mirror fidehties shown in 
figure [2] can be accessed. The spontaneous scattering rate is related to the intensity I{t) 
and detuning A through the expression 

where F is the natural linewidth of the transition, /sat is the saturation intensity and 
I(t) = /exp(— t^/2r^), since the time-dependence of the intensity is exactly that of the 
two-photon Rabi frequency. Note that the Bragg mirror is created from two counter- 
propagating lasers of intensity /, and hence the spontaneous scattering rate has an 
effective intensity of 21. The detuning can be written in terms of the two-photon Rabi 
frequency and intensity via the relations 

_ {dE/hf _ I 

2A~ 2 A ~n%ocA' ^ ' 

where d is the electric dipole moment of the atom, c the speed of light in vacuum and 
^0 the permittivity of free space. Thus the proportion of atoms lost due to spontaneous 
emission is 



si^i.Tj) = / R{i{t),n,T)dt 

/ — dt (21) 

J l + 2Ie-^V^-'/h^, + A{d^/h?e,cY{P/Tm^) ' ^ ' 



By substituting the optimal pulse parameters (cf. figure |3]) into equation (21), we 
can calculate the loss due to spontaneous emission as a function of intensity for any 
order maximum fidelity tt pulse and source momentum width. As an example, consider 
a state-of-the-art laser system (such as that used in [TT]) that addresses the 780 nm 
D2 transition of a source of '^''Rb atoms. This transition has a natural linewidth of 
r = 27r X 6.07 MHz, an electric dipole moment d = 2 x 10~^^ C/m, and a saturation 
intensity of /sat = 1-6 mW/cm^. These numbers have been used to calculate the curves 
in figure [Sj which shows S{ilopt, Topt, I) for a range of Bragg orders. Note that a loss over 
100% indicates that all the atoms in the cloud have spontaneously scattered a photon 
before the pulse has finished. 

Now suppose that our experiment can tolerate no more than 1% of the atoms 
lost due to spontaneous emission, and furthermore we are limited to a laser intensity 
of 1 W/cm^. Then the shaded region in figure [s] indicates the Bragg orders in our 
example scenario where the maximum mirror efficiency can be achieved. Specifically, 
maximum efficiency Bragg mirrors can certainly be achieved for n < 5; however they 
are unobtainable for n > 10. A maximally efficient 10th order Bragg mirror could be 
achieved by doubling the maximum available laser intensity to 2 mW/cm^. 

Figure |5] illustrates that even with a state-of-the-art apparatus, maximum fidelity 
mirrors are practically unobtainable for moderate Bragg orders. This begs the question: 
how does a clamp on the two-photon Rabi frequency effect the maximum mirror fidelity? 
Towards answering this question, we numerically calculated some maximum fidelities for 
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Laser intensity (units of W/cm^) 



Figure 5. Plots of the percentage of atoms lost due to spontaneous emission as a 
function of the laser intensity for a range of Bragg orders. The curves shown have 
been calculated using the optimal pulse parameters shown in figure [3] and are for a 
a = 0.1 hk cloud; there is little change for other momentum widths less than Ihk. 
The shaded rectangle indicates at what Bragg orders the maximum fidelity can be 
achieved, assuming that the interferometer cannot tolerate greater than 1% loss due 
to spontaneous emission, and that the maximum laser intensity is 1 W/cm^. 
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Figure 6. Plots of the maximum fidelity as a function of momentum width for 
a range of Bragg orders. The optimum two-photon Rabi frequency was subject to 
the constraint flopt < i^max, where fimax = 20 Ur- For n > 4, the optimum il that 
gave the maximum fidelity was always ilopt = i^max- As before, the numerical data is 
represented by the points, and the lines are simply to guide the eye. 



a two-photon Rabi frequency clamp of fimax = 20 cOr- The resuhs of this investigation are 
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shown in figure [6j At n = 3, the optimal two-photon Rabi frequency Qopt is always less 
than ^max; and so the clamp has no effect on the maximum fidelity. However, for n > 4, 
^max is less than the optimal two-photon Rabi frequency. This constrains the allowable 
optimum to f2opt < ^max, the effect of which is to drastically reduce the fidelity. In 
particular, figure [6] illustrates two points. Firstly, when the two-photon Rabi frequency 
is constrained increases in n result in dramatic decreases in the maximum obtainable 
fidelity. The second point to note is that larger Bragg orders require an increasingly 
narrow momentum width in order to obtain good transfer efficiencies. For example, at 
n = 3 a fidelity greater than 0.9 can be achieved for any a < 0.1 hk. For n = 5 and 
n = 6, the momentum width needs to be a < 0.015 hk and cr < 0.001 hk, respectively, 
in order to achieve a Bragg mirror with fidelity greater than 0.9. 

Clamping the the two-photon Rabi frequency changes both the optimal two-photon 
Rabi frequency and the optimal pulse duration. When fimax = 20 Ur, the optimal pulse 
time increases from Topt ~ 0.4 — 0.7 for n = 4 to Topt ~ 2 — 2.7 co~^ for n = 5 and 
Topt ~ 19 — 26 uj~^ for n = Q. For the ^''Rb 780 nm D2 transition, this translates to pulse 
durations on the order of 1-10 ms. For pulse times considerably longer or shorter than 
these optimal values, the fidelity will be close to zero. 

To summarize, limits on laser power and spontaneous loss place an upper bound 
on the two-photon Rabi frequency. Clamping the two-photon Rabi frequency reduces 
the efficiency of higher order Bragg processes by a factor that strongly depends upon 
the momentum width of the cloud. As n increases, this degradation in fidelity becomes 
more severe and can only be overcome by increasingly narrow momentum widths and 
longer pulse durations. 

We now calculate the maximum fidelities that can practically be achieved for 
the experimental scenario considered when calculating figure |5| For sufficiently large 
detuning, we can ignore the time-dependence in the denominator of equation (21). In 
this case 

~ 47r 1 + 2///,at + 4(rfVn%oc)2(/V^2^]2) • (22) 

Rearranging this expression gives, for a fixed S, the maximum two-photon Rabi 
frequency function of / and r: 

a„„(,,.,.^ja£M^=. (23) 



'{rr/(v^S)- 2} 

Figure [7](a) shows the effect of this constraint on the maximum fidelities for S = 0.01 
and a cloud of momentum width a = 0.1 hk. It is clear from this plot that a maximally 
efficient Bragg refiection with F.,^ ^ 0.92 (cf. figure |2]) requires some minimum laser 
intensity /min, and that this minimum intensity increases with increasing Bragg order. 
For intensities less than /min, the efficiency of the Bragg mirror is vastly reduced. 
Figures [7](b) and (c) show the optimal pulse durations and detunings, respectively, 
needed to achieve the constrained maximum fidelities of figure [7|^a). 
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Figure 7. (a) Plots of the maximum fidelity as a function of laser intensity for 
a cr = O.lhk cloud, where the two-photon Rabi frequency has been constrained by 
equation (23 1. For each Bragg order, the optimal pulse duration and detuning are 
given as functions of intensity in (b) and (c) respectively. The lines are simply to guide 
the eye. 



5. Analysis for Mach-Zehnder interferometer 

In section [4], we used loss from the target state as a characterization of Bragg mirror 
efficiency. However, the phase of the atoms also has an important effect on the SNR. 
We therefore consider the effect of losses and interference by extending our analysis to 
a MZ interferometer. We show that the MZ interferometer shares many of the trends 
exhibited by the Bragg mirror, including limitations on the SNR due to the momentum 
width. Furthermore, we calculate the optimal Bragg order and (for a thermal source) 
momentum width that maximize the SNR. Finally, we compare the SNR for optimized 
thermal and Bose-condensed sources. 



5.1. Efficiency criterion for Mach-Zehnder interferometer 

For a MZ interferometer, we ultimately want to maximize the SNR Q for a phase 
measurement. A good figure of merit will be proportional to the SNR, and will 
incorporate losses and phase information. Losses are accounted for by simply counting 
the number of atoms in the two output states of the interferometer. Phase information 
is related to the SNR via the fringe contrast of the relative populations. We thus define 
the MZ interferometer efficiency criterion as 

G = Cy/N. (24) 

It is not difficult to generalize the contrast C and the atomic population N in the two 
output ports of the interferometer to the case of finite momentum width atomic clouds. 
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Define the total population of momentum state mhk as 



P^{<P) = [ -L=e-'''l^^'' |c^(5,0,t/)|' db, (25) 



where is now the time after all pulses have interacted with the atomic cloud. Without 
loss of generality, we normalize the total number of atoms to 1. In this case, the 
coefficients Cm(5, 0, t) are calculated from the set of equations ([s]) with pulse sequence 
([T]). The relative total population in state nhk is then 



The contrast is given as the difference between maximum and minimum relative total 
populations: 

C= max P'^\4>)- min P'^\4>). (27) 

0<</><27r/n 0<</.<27r/n 

Note that for an nth order Bragg process, the phase is imprinted on the atoms n 
times. Hence the period of any fringes will be 2TT/n. If the fringes are approximately 
sinusoidal, which appears to be true even for moderately large momentum widths, then 

maxo<<^<2./nPr'(0) = Pn\T^/n) and mino<<^<2./n ^^'(0) = ^^'(0)- We assume the MZ 
interferometer operates at mid-fringe, and it is at this point that we determine the total 
number in the two output ports: = P_„(7r/2n) + P„(7r/2n). Therefore the G-factor 
is calculated via the expression 

G ^ ( max Pf{<P) - min Pl^\<P)] ^PlJ^A- (28) 

\ 0<<p<27r/n \}<<p<2iv/n I 

5.2. Two-level model for the Mach-Zehnder interferometer 



As in section 4.2, we can gain qualitative understanding by considering the behaviour 



of an MZ interferometer composed of square pulses operating in the Bragg regime. The 



solution to (14a) and (146) can be expressed as the unitary matrix 



f/(t,0) 



cos (^t) - sin (^i) sin 

^ sin (^t) cos (^t) + sin 



where we use the dimensionless units x = 4urn6 /Qcs and i = Qcst- The phase 
factor exp(±m0) is due to the laser pulse phase, 0. Formally, this is included by 



the replacement f2eff — ^ f^efrexp (m0) in (14a) and (146). For the above unitary matrix 



7r/2 and vr pulses are given by setting t = n /2 and t = vr, respectively. 

During the time between pulses, the atoms undergo free evolution as described by 

We can see then that the output of a MZ interferometer with three pulses, with an 
overall phase shift included on the third pulse, is given by 

iV^output) = f/(7r/2,0)f/^(T)[/(7r,O)f//(r)f/(7r/2,O)|^i„p,t). (31) 
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Figure 8. Plot of the G- factor as a function of w. The key feature to note is that G 
decreases with w (i.e. decreases with increasing momentum width a or Bragg order n). 
However, one can decrease w to obtain the same G-factor by increasing the effective 
Rabi frequency. 



Here T is the interrogation time between pulses. We choose the initial state iV'input) such 
that all the atoms are initially in the ground state, l^/'output) can be used to find the total 



population of the —nhk and nhk states via (25). These can be used to calculate the 



G-factor (28). Provided T is much larger than the pulse duration, G will only depend 



on the w parameter defined in (18). A plot of G as a function of w is shown in figure [8j 
The same conclusions that applied to the Bragg mirror also apply to G, and hence the 
SNR of the interferometer: (a) larger momentum widths and/or higher Bragg orders 
decrease G, and (b) decreases in G due to increases in Bragg order can be neutralised by 
increasing the effective Rabi frequency. However the same limitations on the two-level 
model of the Bragg mirror also apply to this model for the MZ interferometer: it is only 
vahd in the regime a <^hk and ficfr ^ Suriji — l)"/[(ra — 1)!]^. 

5. 3. Numerical simulation for a Mach-Zehnder interferometer 

We performed a numerical optimization for a pulse sequence in the MZ configuration. 
This problem requires the optimization of beamsplitter and mirror pulses, parametrized 



by (f2bs,Tbs) and {^m,Tm) respectively, such that the G-factor (28) is maximized. The 
optimal amplitude and duration of the mirror pulse was chosen based on the values 
calculated by the independent optimization of the Bragg mirror considered in section |4] 
The interrogation time T was always chosen to be much larger than any of the pulse 
durations. In this regime, we expect G to be approximately independent of T. Thus 
the optimization of the MZ interferometer reduces to an optimization of the amplitude 
and duration of the beamsplitter pulses. 

The numerical results for the optimization of G as a function of n and cr with 
unconstrained Q are shown in figure |9| As with the Bragg mirror, it is clear that 
increasing the momentum width leads to a reduction in the maximum G obtainable. 
The other interesting feature is that, apart from the sharp drop in maximum G for the 
lowest two or three n, the maximum G-factor is approximately constant with increasing 
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Figure 9. Plots of the G-factor as a function of Bragg order, n, for different 
momentum widths. RecaU that the total atom number has been normalized to 1. 
No constraint has been placed on the two-photon Rabi frequency. The error bars are 
estimates of the uncertainty in the maximum G-values. Lines joining points are just 
to guide the eye. 



Bragg order, at least for the Bragg orders and momentum widths shown here. Again 
the MZ interferometer behaves similarly to the Bragg mirror; with ample laser power 
the interferometer SNR does not appreciably decrease with increasing Bragg order. 

Whilst G should be independent of the interrogation time T in the limit T — )■ oo, 
computational limitations sometimes necessitated that the simulations be run with 
an interrogation time somewhat smaller than the interrogation times used in typical 
experiments (although T ^ t^, Tbs was always satisfied). This meant that our numerical 
calculations of maximum G had a slight T dependence. Calculations of maximum G for 
different values of T allowed us to estimate the maximum G in the limit T — > oo, and 
place uncertainties on this estimate. The points and error bars in figure|9]are an estimate 
of the maximum G-factor in the limit T — t- oo and the uncertainty on that estimate, 
respectively. Maximum G and associated uncertainties were calculated similarly for the 
remaining plots in this paper. 

We noted in section the approximate quadratic scaling of the optimal mirror 
pulse amplitude with Bragg order. Similarly, the optimal beamsplitter pulse scales 
approximately quadratically with n. Furthermore, the optimized values of (f2bs,Tbs) 
did, in general, differ to the optimal values of {Qm, Tm), but by no more than a factor of 
two (and usually much less). 

We also considered the effect of a two-photon Rabi frequency clamp on the G- 
factor. Clamping at ^rnux = 20 a;^ resulted in significant decreases in the G-factor 
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Figure 10. Plots of riGmax as a function of Bragg order, n, for different momentum 
widths. The two-photon Rabi frequency has been clamped at ilmax — 20 w,.. Note that 
for each a there exists an optimal Bragg order n that maximizes nGmax (and therefore 
the SNR). Lines joining points are just to guide the eye. 



for n > A. Whilst such decreases were amehorated by a narrower momentum width, 
overall the maximum G decreased for increasing n. More interestingly, we compared 
how a clamped two-photon Rabi frequency affected nGmax (which is proportional to the 
SNR), where Gmax is the maximum G-factor for a particular a and n. As shown in 



figure [10| for a fixed a there exists a Bragg order that maximizes nG^ax- Furthermore, 
this optimum Bragg order appears to decrease as a increases. This is a key result, for in 
practice there will always exist some upper bound on the two-photon Rabi frequency. 
Consequently, it is simply not true that increasing the momentum transferred by the 
beamsplitters and mirrors will always improve the SNR, even at the shot noise limit. 

5.4- Comparison of thermal and Bose-condensed sources 



Naively, one may think that figure [TO] implies that a narrower momentum width gives a 
larger SNR for a MZ interferometer. This is only true if the atom flux is fixed for different 
momentum widths. In reality, broader sources tend to have a higher flux than narrower 
sources. One therefore expects there to be some trade-off between increasing the G- 
factor by narrowing the momentum width, and increasing the atom flux by broadening 
the momentum width. In this final section, we include the flux considerations due to 
choice of source in order to compare the SNR of thermal and Bose-condensed sources. 
We can incorporate the Bragg order and the effect of source type on flux into our 



Why momentum width matters for atom interferometry 



21 



analysis by defining the effective G-factor: 

Geff = n^/W^G^^,, (32) 

where Ni is the flux of the source that is input into the interferometer. Ges is directly 
proportional to the SNR, and includes considerations due to the initial flux of the source, 
momentum transferred by the beamsplitters and mirrors, phase information (via the 
fringe contrast) and losses from the two interferometer modes during the interferometric 
process. 

Thermal sources are typically velocity-selected before the flrst beamsplitting pulse. 
Assuming a velocity-selection pulse reduces the cloud to a Gaussian distribution of width 
a, then the fraction of atoms remaining after the velocity-selection pulse is proportional 
to cr/ao, where (Tq is the initial momentum width of the thermal source. We therefore 
assume Ni = a/ao for a thermal source, where for convenience we have normalized the 
flux of the initial thermal source of width ctq to 1. By maximizing Gcs we can not only 
determine the optimal Bragg order, but also the optimal velocity-selection. 

Unhke a velocity-selected thermal source, the input flux of Bose-condensed sources 
does not linearly depend upon the momentum width. BEG and atom laser fluxes 
typically depend upon a multitude of factors that are highly apparatus dependent. For 
the purposes of comparing Ges for thermal and Bose-condensed sources, we assume the 
flux of an expanded BEG or atom laser to be a factor of 25 lower than a Ihk thermal 
source. This is the difference in flux between the state-of-the-art BEG machine in [T7] 
and the best published state- and velocity-selected thermal cloud used in a gravimeter 
|18j . We also assume an expanded BEG and atom laser have momentum widths of 
0.1 hk [12] and 0.01 hk [27], respectively. 

Figure 11 graphically compares the effective G-factor for thermal and Bose- 
condensed sources. For a thermal source, these data indicate that a better SNR is 
obtained by velocity-selecting less (giving a larger initial flux) and operating at a lower 
Bragg order where SNR reductions due to the broader momentum width and limited 
laser power are reduced. Indeed, the optimal pulse amphtudes and times for n = 3 are 
unaffected by the f^max = 20 ojr clamp we have considered, whilst the clamp drastically 
reduces Gmax for the 5th and 6th Bragg orders. In contrast, the narrow momentum 
widths of Bose-condensed sources allow for low loss, high contrast operation at higher 
Bragg orders than any thermal source. For the state-of-the-art parameters consider 
here, the best SNR for an expanded BEG and atom laser occur at n = 4 and n = 5, 
respectively. The theoretical SNR of a plane-wave Bose-condensed source, also shown in 
flgure 11, indicates that further narrowing of the momentum width of Bose-condensed 
sources allows for operation at increasingly larger Bragg orders. It appears, therefore, 
that LMT atom-optical elements are more beneflcial for increasing the SNR for Bose- 
condensed sources than for thermal sources. 



From figure 11 , we can see that current state-of-the art thermal and Bose-condensed 
sources theoretically give, within a factor of two, the same SNR for a MZ interferometer. 
This shows that with the application of LMT beamsplitters and mirrors, Bose-condensed 
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Figure 11. Plots comparing the effective G-factor for thermal and BEC sources. 
The plot on the left-hand side displays, for different Bragg orders, Gee as a function of 
momentum width for a thermal source initially of width = Ihk velocity-selected to 
momentum width a. The plot on the right-hand side shows Gofi as a function of Bragg 
order for an expanded BEC source of a = 0.1 hk and an atom laser of cr = 0.01 Kk. The 
dotted line in the right-hand figure is Goff in the limit of a plane- wave Bose-condensed 
source. We have assumed the Bose-condensed sources to have a factor of 25 less flux 
than the Ihk thermal source. The scale for Ges is the same for the thermal and Bose- 
condensed sources, thus allowing for easy comparison. These data assume an upper 
bound on the two-photon Rabi frequency of f2max — 20 cj^. 



sources are competitive with thermal sources, even with a factor of 25 less flux. 
Technological improvements that increase the flux of Bose-condensed sources are not 
inconceivable, and could be achieved by increasing the atom number of the condensate 
and/or increasing the repetition rate of the BEC machine. For the situation considered 
in figure 11, an atom laser with a comparable flux to a, Ihk thermal source would give 
a SNR four times larger than the best SNR achievable with this thermal source. It is 
likely that Bose-condensed sources have further advantages when the momentum width 
transverse to the direction of the optical lattice is considered. Our analysis entirely 
ignored the transverse momentum width, and consequently effects on the SNR due to 
wavefront distortions and the Coriolis effect were neglected. Based on work by Debs et 
al. [12] and Louchet-Chauvet et al. [28], we expect narrower momentum width sources 
to be less susceptible to aberrations and the Coriolis effect. 



6. Conclusions 



We have developed a relatively straightforward theoretical method for optimizing Bragg 
atom-optical elements. We have used this method to investigate the effect of momentum 
width on the efficiency of Bragg mirrors and Bragg MZ atom interferometers. Our 
comprehensive numerical analysis of the Bragg mirror (MZ interferometer), coupled 
with the qualitative understanding given by a two-level model, demonstrated that (a) 
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the source momentum width fundamentally hmits the efficiency of Bragg mirrors (SNR 
of MZ interferometers); (b) for fixed momentum width, increasing the Bragg order does 
not generally change the maximum fidelity (G-factor) provided one has sufficient laser 
power; (c) the optimal two-photon Rabi frequency for mirror and beamsplitter pulses 
scales roughly quadratically with the Bragg order; and (d) a limitation on the two- 
photon Rabi frequency will, above some Bragg order, result in drastic decreases in the 
maximum mirror efficiency (G-f actor), and further restricts a high fidelity (G-f actor) 
to narrow momentum width clouds. Finally, we incorporated source type and, for a 
particular two-photon Rabi frequency clamp, calculated the optimal Bragg order and 
(for a thermal source) velocity selection that maximized the SNR. We compared thermal 
and Bose-condensed sources, and showed that state-of-the-art Bose-condensed sources 
are competitive with state-of-the-art thermal sources. 

The analysis undertaken in this paper shows that, fundamentally, momentum 
width does matter for atom interferometry with Bragg pulses. Moreover, we have 
demonstrated that this analysis can be used to optimize the atomic source, beamsplitters 
and mirrors for a wide range of Bragg diffraction experiments and interferometers. 
For example, one could, in tangent with a few practical constraints, use our analysis 
to calculate the optimal velocity-selection pulse to use for a thermal source MZ 
interferometer. Given the potential to increase BEG fiux, alongside other practical 
advantages of a narrow momentum width such as a lower susceptibility to wavcfront 
distortions and the Coriolis effect, we believe Bosc-condcnscd sources to be a viable 
candidate for the next generation of precision sensors based on Bragg pulse atom 
interferometry. 
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